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Abstract 

The study of the rotational spectrum of NaCN (X'A') has recently been extended in frequency and in quantum numbers. Dif- 
ficulties have been encountered in fitting the transition frequencies within experimental uncertainties. Various trial fits traced the 
difficulties to the incomplete diagonalization of the Hamiltonian. Employing fewer spectroscopic parameters than before, the tran- 
sition frequencies could be reproduced within experimental uncertainties on average. Predictions of a-type R -branch transitions 
with K a < 7 up to 570 GHz should be reliable to better than 1 MHz. In addition, modified spectroscopic parameters have been 
derived for the 13 C isotopic species of NaCN. 
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i . There are several ways to form a molecule from a metal M 
O and the CN group. Several metals, such as the transition met- 
als copper and zinc, prefer a linear cyanide arrangement with 
j-^j an MC bond. Linear isocyanides with an MN bond are pre- 
'— 1 ferred by others such as the main group metals magnesium and 
| aluminum. Sodium and potassium, however, display a different 
structural motif, a nearly T-shaped structure with almost equally 
(**T} long MC and MN bonds. As a consequence, both molecules are 
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asymmetric top rotors. 

The rotational spectrum of NaCN has been studied by van 
Vaals et al. [1] between 9 and 40 GHz employing molecular- 
beam electric -resonance spectroscopy. 20 rotational transitions 
were recorded, half of them being much weaker fo-type tran- 
^ sitions. NaCN was detected in the circumstellar envelope of 
• • the carbon-rich star CW Leo, also known as IRC +10216 J2l, 
. £h soon a f ter th ree pairs of related unidentified lines, also detected 
in that source, had been attributed to MgNC fl. Recently, 
$— i He et al. |4] described a molecular line survey of CW Leo at 
5^ 2 and 1.3 mm. In order to analyze the NaCN emission, the 
spectroscopic parameters of NaCN were reevaluated because 
initially, the quartic centrifugal distortion parameters were ex- 
pressed as r values, which are rarely implemented in current 
spectroscopy programs. The predictions derived from the re- 
sulting parameters were good enough for the NaCN emission 
features obtained in the course of that line survey as well as 
for previously published transitions. Observations at slightly 
higher frequencies and with a better signal-to-noise ratio, how- 
ever, revealed increasing deviations largely with increasing J, 
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prompting Halfen and Ziurys ^ to study the rotational spec- 
trum of NaCN between approximately 180 and 530 GHz. They 
encounterted difficulties in reproducing their experimental tran- 
sition frequencies on average within experimental uncertainties 
and attributed this to the floppy nature of the molecule. 

In the course of creating an updated catalog entry for the 
Cologne Database for Molecular Spectroscopy (CDMS~ 
the previously reported experimental data for NaCN were 
critically evaluated. Using fewer spectroscopic parameters than 
before, it was possible to reproduce both data sets well if the 
Hamiltonian was diagonalized sufficiently. 

2. Considerations for the spectral analysis 

One finds rather frequently that the set of spectroscopic pa- 
rameters, which reproduces a given list of transition frequen- 
cies within experimental uncertainties, is not unique. One way 
of reducing the ambiguity of a parameter set is to reduce the 
rms error, also knows as the weighted rms, of the fit with as few 
parameters as possible and reasonable. A reasonable parame- 
ter set requires certain lower order term to be used before one 
of higher order is allowed to be included; e.g., the inclusion of 
Ljjk in the fit requires that at least one of the two parameters 
Hj or Hjk has already been used in the fit. In the absence of 
strong correlation among the parameters, it suffices to search 
for the one which reduces the rms error the most. However, 
such correlation is rather common, such that searches for small 
parameter sets may be more complex. It is frequently assumed 
that the smallest parameter set, which reproduces the parame- 
ter set best, will provide the most reasonable predictions. It is 
useful in this context to address the issue of uncertainties at- 
tributed to the transition frequencies. Besides appropriate esti- 
mates and, less frequently in recent years, no estimates at all, 
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one also finds uncertainties which seem to be too optimistic or 
too pessimistic, leading to weights in enlarged data sets which 
are too high or too low, respectively, if these estimates stay un- 
corrected. If the list of transition frequencies is sufficiently large 
with respect to the number of spectroscopic parameters varied 
in the fit, an rms error of 1 .0 overall as well as for well defined 
sub-sets of the frequency list can be viewed as an indication of 
reasonable uncertainty estimates, though caution is advised if 
a very small set of the lines depends particularly strongly on 
one or few parameters. An additional criterion for reasonable 
uncertainty estimates is the distribution of the residuals. Most 
of the lines, about two third, should fit within the uncertain- 
ties, but there may be some lines with residuals exceeding three 
times the uncertainties, however, their number should decrease 
very rapidly with the magnitude of the ratio between residual 
and uncertainty. It is possible that the overall uncertainty es- 
timates in a line list are too pessimistic, but one or few lines 
have residuals which are too large. A rejection threshold is ap- 
plied in final fits in particular in larger data sets as a small set 
of transition frequencies with large residuals may complicate 
the determination of higher order parameters. A threshold of 
3<x is rather common, see e.g. Ref. i8j,|9J], but other thresholds, 
usually less strict, may also be found. 

3. Spectral analysis 

As mentioned in the introduction, NaCN as well as KCN are 
asymmetric tops, in contrast to most other molecules formed 
from a metal atom and the CN group. Since for both the N atom 
is closer to the metal atom HIES], tne y should be viewed as iso- 
cyanides rather than cyanides. However, following the conven- 
tion, we will continue calling these molecules cyanides. Ray's 
asymmetry parameter k is -0.9567 for NaCN, the molecule is 
thus fairly close to the prolate symmetric limit. It has a large a- 
dipole moment component of jj. a = 8.85 D and a much smaller 
^-component of ~0.2 D0 [4] . In spite of the small ^-component, 
van Vaals et al. reported 10 fo-type transitions with AJ = +1, 
AK a = +1, and AK C = +1 along with 10 a-type transitions with 
AJ = 0, +1, AK a = 0, and AK C = +1; the quantum numbers 
extended to J and K a of 16 and 3, respectively. The b-type tran- 
sitions contained essentially all of the AK a + information in 
that data set as well as by far the most AJ + information be- 
cause most of the a-type transitions followed AJ — selection 
rules. Inspection of their residuals suggested that the reported 
uncertainties are rather conservative; they were assumed to be 
3cr uncertainties in Ref. [4]. Employing thirteen parameters 
and the adjusted uncertainties, the latter authors were able to 
reproduce these lines with an rms error of 0.840. The fairly 
large number of parameters compared to the number of lines 
make it less clear whether the adjusted uncertainties of the tran- 
sition frequencies were estimated appropriately. Nevertheless, 
the data permitted good to reasonable predictions well into the 
millimeter wave region. 
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Very recently, Halfen and Ziurys published an extensive set 
of NaCN transition frequencies all pertaining to the strong a- 
type /^-branch transitions with 11 < J" < 32 and < K a < 6 
between 179 and 531 GHz Jill. In order to provide an up- 
dated entry for the CDMS [6, 7], these data were subjected 
to combined fits together with the previous data QJJ] - Since 
Halfen and Ziurys encountered difficulties in reproducing tran- 
sitions mainly with higher K a values, it was decided to add the 
new transition frequencies with each K a series separately in the 
present analysis. Moreover, a tight rejection threshold of 3cr 
was applied since it seemed plausible that large residuals of one 
or few lines may have caused the difficulties. Pickett's spcat 
and spfit programs [11] were used for prediction and fitting of 
the spectra, respectively, as has been done in the two previous 
investigations 0]. 

Starting with the parameter set of Ref. 04J], the previous tran- 
sition frequencies plus 12 new ones with K a — could be re- 
produced overall within uncertainties by adding hi and omitting 
Hk and /i3. At this stage, only one of the previous lines deviated 
from the calculated frequency by slightly more than three times 
the uncertainty. In later fits, the amount and identity of the lines 
displaying larger residual changed sometimes. 

The sequential addition of transition frequencies having K a — 
1-5 afforded one new parameter each to be included. These 
were h\, /13, Ljjk, h, and Lkkj- The addition of K a = 6 lines 
did not require any new parameter. Trial fits with Hk included 
in the fit yielded a value of 41 + 55 Hz with changes in the re- 
maining parameters within the uncertainties. The uncertainties 
of the predominantly /f-dependent terms increased markedly, 
but in no case by more than a factor of 2. Since Hk was not de- 
termined with significance, it was omitted from the final fit. It is 
noteworthy that the uncertainty is a factor of 4 smaller than the 
uncertainty derived in Ref. [4] from the microwave data only. 

This parameter set is smaller than the one in Ref. [5], but 
reproduces almost all of the data within the rejection thresh- 
old. Two of the new lines had residuals around 200 kHz, and 
one line deviated from its calculated position by more than 
300 kHz. Most of the other lines seemed to have been judged 
conservatively with an uncertainty of 50 kHz because transi- 
tion frequencies with K a < 3 gave partial rms errors of less 
than 0.5; the values for K a = 4 and 6 were around 0.6, and 
only lines with K a — 5 yielded a partial rms error of 1.19. The 
previous data yl now yielded a partial rms error of 1.17, indi- 
cating that the previously reported uncertainty estimates have 
been corrected slightly too far. But since none of these lines 
showed particularly large residuals with respect to the uncer- 
tainties, none of these lines were omitted and the uncertainties 
were left unchanged. 

The final set of spectroscopic parameters is given in Table Q] 
together with NaCN and KCN values from Ref. [4] and results 
from a fit of the NaCN microwave data to the KCN parameter 
set. The latter fit has 2 fewer parameters than the one published, 
but has a larger rms error of 1.055 compared with 0.840; this 
was the reason to discard this fit and to choose the one which 
required more parameters. 

The list of transition frequencies with assignments, uncer- 
tainties, and residuals between experimental frequencies and 
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Table 1 : Spectroscopic parameters" (MHz) of NaCN (X 1 A' ) from the present combined fit in comparison to fits employing microwave data only and in comparison 
to KCN along with dimensionless rms errors of each fit. 



Parameter 


NaCN combined fit 


Na 13 CN 


NaCN, microwave data only 


KCN 




Present 


Present 


Ref. [41 


Present 


Ref. \M 


A - (B + C)/2 


50101.74547 (110) 


48021.93283 (314) 


50101 74757 M661 


^01 01 74606 nut 


53528.08107 (268) 


(B + OI2 


7820.189927 (145) 


7652.415049 (391) 


7820.190252 (259) 


7820.190351 (206) 


4737.751389 (93) 


(B-Q/4 


274.1524216 (124) 


272.7985489 (396) 


274.1525055 (177) 


274.1525120 (172) 


100.959727 (64) 


D K X 10 3 


194.11 (25) 


268.89(142) 


196.52 (169) 


193.56 (69) 


956.79 (247) 


D JK X 10 3 


800.06 (2) 


738.90 (39) 


800.23 (34) 


800.24 (19) 


374.609 (50) 


Dj X 1() 3 


13.3123 (5) 


13.2279 (180) 


13.3361 (160) 


13.3313 (96) 


5.24462 (209) 


dy X 10 3 


-2.34537 (23) 


-2.42239 (112) 






-0.51452 (78) 


d 2 x 10 3 


-1.407340 (56) 


-1.379543 (295) 


-1.408951 (182) 


-1.405461 (226) 


-0 237559 ("1991 


H K x 10 3 






1.428 (220) 






H K j x 10 6 


327.4(18) 


303.5 (55) 


155.9(184) 


314.5 (28) 


101.2(11) 


H JK x 10 6 


-26.31 (8) 


-23.43 (30) 


-3.54(120) 


-26.72 (19) 


-7.28(11) 


Hj x 10 9 


-3.77 (34) 


-3.5* 


0.54 (5) 






in x io 9 


-35.68 (17) 


-34.0* 








h 2 x 10 9 


-25.40 (9) 


-24.6* 




-47.65 (199) 


-5.46 (74) 


h 3 x 10 9 


19.34 (9) 


19.1* 


39.92 (688) 






Lkkj x 10 9 


-0.384 (40) 


-0.331* 








Ljjk x 10 9 


-0.845 (43) 


-0.759* 








U x 10 12 


-0.587 (80) 


-0.575* 








rms eiTor 


0.746 


0.807 


0.840 


1.055 


0.827 



"Numbers in parentheses are one standard deviation in units of the least significant figures. 

*Kept fixed to values from the combined NaCN fit multiplied with ratios of the appropriate powers of A — (B + C)/2, B + C, and B — C, see sectionfj] 



those calculated from the final set of spectroscopic parameters 
is given in the Supplementary material. Predictions of the ro- 
tational spectrum along with a documentation will be provided 
in the catalog sectioro of the CDMS. Line, parameter, and fit 
files, along with further auxiliary files, will be available in the 
archive sectior0. 

Ref. |4j also provided spectroscopic parameters for Na 13 CN. 
Since predictions for the main isotopic species from that study 
turned out to be unsatifactory in the upper millimeter and lower 
submillimeter wave regions, it is rather likely that this applies 
to the predictions for Na 13 CN also. 17 rotational transitions 
had been studied by van Vaals et al. |Q]]. The quoted uncertain- 
ties for the derived hyperfine free transition frequencies were 
assumed to be 3 cr values, as has been done for the main iso- 
topolog. Centrifugal distortion parameters were estimated by 
multiplying the NaCN values from the combined fit with ratios 
of the appropriate powers of A - (B + C)/2, B + C, and B - C; 
e. g., H KJ with the ratio of (A-(B + C)/2) 2 (B + C) and h 2 with 
the ratio of (B + C)(B - C) 1 . While this scaling can not be ex- 
pected to hold strictly, it is usually a reasonable asssumption, 
see e. g. Ref. for isotopic species of SiC2 and for Na 13 CN 
or Ref. Hill for the example of the isotopic species of vinyl 
cyanide. 

An rms error of 1 .5 1 8 was achieved by just varying the rota- 
tional and quartic centrifugal distortion parameters. Significant 
reduction of the rms error below values of 1 .0 were obtained 
when either Hkj, Hjk, Hj, or were varied in addition, how- 
ever, the changes for the latter two appeared to be too large 
in comparison to the respective values and were discarded. A 
slightly better fit with an rms error of 0.807 was obtained when 
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both Hkj and Hjk were released in the fit in addition to the 
lower order parameters. The resulting fit is also given in Ta- 
bleQ] Updated predictions forNa 13 CN as well as auxiliary files 
will be provided in the CDMS. 

4. Discussion 

It may be surprising that in the present study a better fit was 
obtained with 17 parameters which is less than 23 employed 
previously J5|]. The reason turned out to be rather mundane. 
The programs spfit and spcat require minimum and maximum 
values for K to be specified. In a prolate type, asymmetric top 
molecule, the minimum value is obviously 0, and the maximum 
value has to be at least as large as the highest K a value in the 
line list. The numerical diagonalization of the Hamiltonian, 
however, usually requires the maximum value to be somewhat 
larger than that value. It appears as if an increase by the high- 
est K a order of the parameters in the fit is a reasonable estimate 
for this increase. Transitions in the NaCN line list extend to 
K a = 6, and the parameter I4 connects levels differing in K a 
by 8, e.g. the K a — 4 asymmetry components. Hence, the 
maximum K value requested in the fit should be around 14. In 
fact, no change in parameters and in the rms error occurs for 
K a > 15. The data set for SO2 in its ground vibrational state 
extends to K a = 28 [18]; trial fits have shown that a maximum K 
value of at least 33 should be requested. Even larger maximum 
K values may have to be requested for predicting a spectrum 
at room temperature for convergence of the calculated partition 
function. 

As can be seen in Table Q] the lower order spectroscopic pa- 
rameters, up to the quartic centrifugal distortion terms, agree 
very well with the previous parameter determined from mi- 
crowave data only. Since merely selected parameters of higher 
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order were employed, in particular for the MW fit, it is prob- 
ably not surprising, that the agreement for the higher order 
parameters is not so good. In the alternative MW data fit, 
which employs exactly the parameters which were used in 
the KCN fit, the parameters Hkj and Hjk agree well with 
the present values, and only /12 differs by almost a factor of 
two. It is interesting to note, that this fit yields much bet- 
ter predictions of the new data [5] than the one published 
earlier [4]. The 19o,i9 _ I80.18 transition, e.g., predicted by 
the alternative MW fit at 284160.864 (461) MHz was mea- 
sured at 284161.701 (50) MHz, almost within twice the pre- 
dicted uncertainty. The previously published fit |01 predicts 
284169.538 (634) MHz. Apparently, the parameter set which 
yields the smallest rms error is not always the best one, and this 
situation is more likely to occur if the number of parameters 
is comparatively large with respect to the number of transition 
frequencies or if a certain parameter is determined by only few 
transitions. 

The present parameter set not only reproduces the experi- 
mental transitions well, but it should also permit to a certain 
degree reliable predictions beyond the experimentally accessed 
quantum numbers. It should be pointed out that the discus- 
sion of the predictive power of the present parameter set is con- 
cerned mainly with a-type transitions because the a-component 
of 8.85 D is very much larger than the ^-component of ~0.2 D. 
The Z?-type transitions are weaker than the a-type transitions by 
(yu a //i/,) 2 « 1960, making them unimportant for radioastron- 
omy. No b-type transitions have been provided in the previous 
CDMS catalog entry, and this will also apply for the new entry. 

We consider a prediction to be reliable if the frequency of 
a transition not yet measured will be found to deviate from 
the predicted frequency by an amount generally not exceeding 
three times the predicted uncertainty and never exceeding that 
value by far. Assuming that the effects of centrifugal distor- 
tion parameters not included in the fit do not exceed the pre- 
dicted uncertainties, predictions may be reliable as long as the 
uncertainties do not exceed 1 MHz. At low values of K a this 
corresponds to / up to around 50. Extrapolation in K a is very 
/-dependent. Around J = 15, K a = 10 may be predicted rea- 
sonably, around / = 35, K a — 9 is the limit, and beyond / of 
40, even predictions with K a - % should be viewed with cau- 
tion. Since the assumption mentioned above is not necessarily 
justified, all extrapolations should be viewed with additional 
caution. 

The parameter Hk may be the most important one not yet in- 
cluded in the fit which will have effects on the calculated tran- 
sition frequencies and energies. Using a value of Hk = 200 Hz 
as fixed, a value almost four times the uncertainty of Hk from 
trial fits, see section [3] and thus near the upper end of pos- 
sible values, changes the energies by less than 1 cirT 1 for 
/ = K a < 21, The a-type transitions frequencies change by 
insignificant amounts for the / and K a range for which the pre- 
dictions are deemed to be reliable, see above. 

The easiest way to improve the Zf-level structure of NaCN 
will be recording higher-/ a-type /^-branch transitions with 
AK a = because of their favorable intensities. The K c = / tran- 
sition are already approaching oblate pairing as the / = 33 - 32 



transition frequencies differ by only 60 MHz. Transitions ap- 
proaching oblate pairing contribute to the determination of 
purely /f-dependent parameters. At still higher /, these tran- 
sitions will reach oblate pairing and transitions with K c — J - 1 
etc. will approach and eventually reach oblate pairing also. In 
fact, the uncertainty of Hk in present trials fits was reduced to 
55 Hz (see section [3), a factor of 4 smaller than the value of 
220 Hz obtained from microwave data only fl, see also Ta- 
ble Q] indicating that the newly measred a-type R -branch tran- 
sitions [5] have improved the purely -dependent parameters. 

Obviously, AK a + transitions probe the A'-level structure 
more directly. In that case, it may be more promising to search 
for AK a = 2 transitions which are frequently much stronger 
than fr-type transitions, in particular in the millimeter wave re- 
gion. Transitions with K a < 6 and not too high / should be 
predicted sufficiently well with the largest deviations from the 
predicted frequencies of order of 10 MHz. 

The centrifugal distortion parameters which have been varied 
in the Na I3 CN fit are close to the initial values obtained by scal- 
ing NaCN values. The most pronounced deviation occurs for 
Dk, which shows a very large isotopic shift. Such large shift, 
however, can be attributed to the accidentially small value of 
Dk compared with D jk and is not uncommon. A large isotopic 
shift was also observed, e. g., for Ak of 37 C1 35 CK)2 compared 
to 35 C1 CIO2 which was fully accounted for by the harmonic 



force field calculation 111 311 . A relatively much smaller deviation 
occurs for d\, whose magnitude for the I3 C isotopoic species is 
actually slightly larger than that of the main species. It is pos- 
sible that this unexpected relation is caused by the assumptions 
made for some of the fixed higher order parameters. 

Predictions of the rotational transition frequencies of 
Na 13 CN should be viewed with much greater caution than those 
for the main isotopic species because of the much smaller data 
set. Nevertheless, the predictions should be reliable to reason- 
able at least up to / of 15, maybe 20, and up to K a of 4, maybe 
5. 

Halfen and Ziurys derived a rotational temperature of 
about 41 K for NaCN in the circumstellar envelope of CW Leo. 
The Boltzmann peak at this temperature occurs slightly above 
200 GHz, and the line intensities appear to be negligible above 
500 GHz. Moreover, transitions with K a - 10 are around 3 
orders of magnitude weaker than those with K a — or 1 . The 
predictions of the NaCN rotational spectrum derived from the 
present study should thus be adequate for all observational pur- 
poses, at least as far as the ground vibrational state is concerned. 
The bending mode (i>2 = 1) is the excited vibrational state low- 
est in energy. It is -170 cirT 1 or ~245 K above ground [Q]]. 
Transitions pertaining to this state will be hard to detect in the 
absence of infrared pumping, which, however, cannot be ruled 
out. 

5. Conclusions 

The rotational spectrum of NaCN has been reanalyzed. The 
previous difficulties in reproducing the experimental transition 
frequencies U have been attributed to the incomplete diago- 
nalization of the Hamiltonian. The predictions generated from 
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the present data set should be good enough for all observational 
purposes restricted to the ground vibrational state. The modi- 
fied parameter set for Na 13 CN should provide better prediction 
of the rotational spectrum up to moderate quantum numbers. 
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Appendix A. Supplementary material 

Supplementary data for this article are available on Sci- 
enceDirect (www.sciencedirect.com) and as part of the 
Ohio State University Molecular Spectroscopy Archives 
(http://library.osu.edu/sites/msa/jmsa_hp.htm). Supplementary 
data associated with this article can be found, in the online ver- 
sion, at doi: . 
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